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We consider dynamically generated spin squeezing in interacting bimodal condensates. We show 
that particle losses and non-zero temperature effects in a multimode theory completely change 
the scaling of the best squeezing for large atom numbers. We present the new scalings and we give 
approximate analytical expressions for the squeezing in the thermodynamic limit. Besides reviewing 
our recent theoretical results, we give here a simple physical picture of how decoherence acts to limit 
the squeezing. We show in particular that under certain conditions the decoherence due to losses 
and non-zero temperature acts as a simple dephasing. 
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I. INTRODUCTION 

Spin squeezing is about creating quantum correlations 
in a many-body system that can be useful for metrol- 



ogy. An example is that of atomic clocks, where the 
ultimate signal to noise ratio, once all the technical noise 
has been eliminated, can be improved by manipulating 
and controlling the system at the level of its quantum 
fluctuations. 



A. Spin squeezing and atomic clocks 

The aim of an atomic clock is to measure precisely the 
energy difference between two atomic states a and h that 
are for example two hyperfine states of an alkali atom. 
To explain how the clock works and to introduce spin 
squeezing, we shall describe the ensemble of N atoms 
used in the clock using the picture of a "collective spin" 
that evolves on the so-called Bloch sphere. The collective 
spin is simply the sum of the effective spins 1/2 that 
describe the internal degrees of freedom of each atom. In 
the second quantized formalism the three hermitian spin 
components Sx, Sy and Sz are defined by: 



Sx + iSy — a^b , 



Sz 



b^b 



(1) 



lOl where and b^ are creation operators of particles in the 
internal states a and b. The spin operators are dimension- 
less and obey the commutation relations [Sx,Sy] — iSz 
and cyclic permutations. For the moment we do not care 
about the external degrees of freedom of the atoms. The 
component Sz of the collective spin is half the popula- 
tion difference between states a and 6, while Sx and Sy 
describe the coherence between these states. If the N 
atoms are prepared in a coherent superposition of states 
a and b with relative phase 20: 



1 



)N = 



V2 



N 



|0) 
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where |0) is the vacuum, the collective spin lies on the 
equatorial plane of the Bloch sphere, pointing at an an- 
gle —20 with respect to the x axis. For non-interacting 
atoms, the further evolution is ruled by the Hamiltonian 



Hq — tkJabSz 



(3) 
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FIG. 1: (Color online) Uncorrelated state and squeezed state 
represented on the Bloch sphere. 



and the spin precesses around the z axis with the Larmor 
frequency ujab- The atomic clock measures the phase ac- 
cumulated by the collective spin during a long precession 
time r. From this phase the frequency tUab is deduced. 
Atomic clocks are nowadays so precise that they are sen- 
sitive to the quantum noise of the collective spin. If for 
example the spin is initially prepared along the x axis (in 
an eigenstate of S'^;), it has necessarily fluctuations in the 
transverse components Sy and Sz such that: 



N 



(4) 



In particular fluctuations of Sy introduce a statistical 
variance on the accumulated phase during the precession 
time and on the measured frequency LUat- In the uncorre- 
lated state ([2]) with = 0, root mean square fluctuations 
of Sy and Sz are equal: ASy — ASz — VN /2. In a Ram- 
sey measurement with interrogation time r, these quan- 
tum fluctuations introduce the root mean square fluctu- 
ations of the measured frequency equal to jl| : 



(5) 



Nt 



This noise coming from quantum fluctuations, intrinsic to 
the initial state where each atom is in a superposition of 
a and b, is known in clocks as the "partition noise" . The 
idea of spin squeezing 0] is that the Heisenberg relation 
(|4]) allows to reduce ASy provided that ASz is increased. 
This idea is illustrated in FiglTJ To quantify the spin 
squeezing we use the parameter introduced in [1] : 



(6) 



where N is the total atom number, AS'j_ is the min- 
imal variance of the spin orthogonally to its mean value 
(S). The state is squeezed if and only if < 1. As 
explained in [l| , ^ directly gives the reduction of the sta- 
tistical fluctuations of the measured frequency uJat with 
respect to uncorrelated atoms, for the same atom number 
N and the same Ramsey time t: 



Nt 



(7) 



The parameter ^ in Eq.(l6]) is in fact the properly nor- 
malized ratio between the "noise" AS'j.^min and the "sig- 
nal" |(S)|. In experiments AS'_L^niin is directly measured 
by measuring Sz after an appropriate state rotation and 
|(S)| is separately deduced from the Ramsey fringes con- 
trast. 



B. State of the art 

On one hand the most precise atomic clocks using mi- 
crowave transitions in cold alkali atoms have already 
reached the quantum partition noise limit with atom 
numbers up to iV = 6 x 10^ 0|. On the other hand, 
very recently a significant amount of spin squeezing, up 
to —8 dB (^^ = 10~° ®) [3] was measured in dedicated, 
proof-of-principle experiments. In Q squeezing was cre- 
ated in a large sample of iV = 5 x 10* atoms with a 
feedback mechanism in a resonant optical cavity, while 
in 4] and in 6] the squeezing was created in smaller 
samples, of order N — 10^, using atomic interactions in 
bimodal condensates. The ultimate limits of the different 
paths to spin squeezing are still an open question. Here 
we concentrate on a dynamical scheme using interactions 
in bimodal condensates 0, [^-Q and analyze in particu- 
lar the influence of dephasing, decoherence and non-zero 
temperature on this squeezing scheme. 



C. Two-mode scalings without decoherence 

We consider for simplicity a bimodal condensate with 
identical interactions in the components a and b with 
coupling constants gaa = gtb = 9 and no crossed a-b in- 
teractions [1^. We assume that the initial state is the 
factorized state ([2]) with (j) = Q and a fixed total num- 
ber of atoms N . In a two-mode picture, interactions 
introduce a Hamiltonian that is non-linear in the spin 
operator: 



F„i = hxSl 



(8) 



The quadratic form ([8]) is obtained expanding the system 
Hamiltonian to second order around the average num- 
bers of particles in components a and 5, Na and iVf,, 
both equal to N/2 for the initial state ^ [i, [I^l- S-X 
is thus the derivative of the chemical potential with re- 
spect to the particle number in each component hx — 
dfia/dNa = d^ib/dNb evaluated in Na — Na, Nb = Nb- 
The general expression of the expanded Hamiltonian in- 
cluding drift terms for non-symmetric interactions, non- 
symmetric splitting or fluctuations in the total particle 
number can be found in [l^, [ll[ . We are interested in 
the best squeezing that can be obtained in the thermo- 
dynamic limit for a spatially homogeneous system: 
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(9) 
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FIG. 2: (Color online) Left: Deformation of quantum fluc- 
tuations due to the non-linear Hamiltonian ([8]). Right: Best 
squeezing time. 



we therefore explicitly write x terms of the interaction 
constant g and the volume V of the system: 



Ma 



gNg 
V 



fj-b 



V 



X 
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(10) 



We can consider the non-linear Hamiltonian (jS]) as a 
Hamiltonian of the form ^ with a Larmor frequency 
iijab that depends itself on Sz- As explained in [3| and 
shown in Figl2] (left), H^i "twists" the transverse spin 
fluctuations and generates spin squeezing. However, in 
order to minimize the evolution should not go to far. 
E.g. when the fluctuations become too much distorted 
and start to wrap around the Bloch sphere, the "signal" 
I (S) I in the denominator of decreases and the squeezing 
parameter increases again. A sketch of the time depen- 
dence of when the state ^ evolves under the influ- 
ence of Hn\ ^ is given in Figl2] (right). We name "best 
squeezing time" t^-m the time that minimizes ^^(t) and 
"best squeezing" ^^j^^ the corresponding squeezing. The 
analytical expression for the squeezing as a function of 
time is given in Q. Introducing an appropriate rescal- 
ing of the time variable as in [l2| , this gives the following 
scalings of the best squeezing and the best squeezing time 
for TV ^ 1 and p,g constants: 



(Cmin)Hni 



3^/3 1 

2 iV2/3 



pgtn 



H„i 



31/6^1/3 ^ 
(11) 

This differs from the original prediction in by numer- 
ical factors. Subsequent studies 0, [Si gave indications 
that the achievement of large squeezing in condensates 
should be possible even in presence of decoherence, but 
were not able to confirm or disprove the scalings ([TT|). 



D. New scalings in presence of decoherence 

We will show in sections lllll and IIVI that the scalings 
(jlip are disproved when decoherence coming from par- 
ticle losses or non-zero temperature is included in the 
description. To summarize, instead of tending to zero 
when N ^ 00 in the thermodynamic limit, S.'^^^ tends to 
a positive constant 



2 

min 



lim. therm. 



constant > 



(12) 



Concerning the best squeezing time imin we distinguish 
two cases. In the case of particle losses tmin is finite in 
the thermodynamic limit and scales as 



oc 



2 

min 



(13) 



In the case of finite temperature tmin cannot be calculated 
within our analytical treatment that neglects interactions 
among Bogoliubov modes. Nevertheless, we introduce 
a "close-to-best" squeezing time t,, defined in equation 
(|37l) . at which approaches £_f^^^ with a finite precision 
T], that scales as 



pgtr, 



1 



oc 



2 

min 



(14) 



provided that trj remains smaller than the typical collision 
time among Bogoliubov modes. 

We will show that these new scalings in presence of 
decoherence are quite general and that the physics of how 
decoherence acts is caught by a very simple dephasing 
model that we shall solve exactly and study in detail in 
the next section. 



II. DEPHASING MODEL 

In this section we consider a dephasing Hamiltonian 
model of the form 



(15) 



where Z? is a Gaussian real random variable of zero mean. 
We assume here that D is time independent, but it varies 
randomly from one experimental realization to the other 
mimicking a stationary random dephasing environment. 
We also assume that D has a variance of the order of N 
for N large H: 



N 



N ■ 



(16) 



Finally enoisci finite in the thermodynamic limit, is a 
small parameter of the theory and we limit ourselves in 
general to first order in this quantity. An exception is 
made in subsection IIIBI where expressions to all orders 
in enoisc are given. 

Starting with the initial state ^ with (f> — we will 
show that this minimal model reproduces the scalings 
and pi)) . In the subsequent sections Hill and HVl we 
will detail an analogy between the dephasing model and 
microscopic models accounting for the effect of particle 
losses or of non-zero temperature on squeezing. The pa- 



to the lost fraction of particles or the populations of ther- 
mally excited modes, respectively. 
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A. Squeezing in the thermodynamic Hmit 

For the symmetric case we consider, the mean spin is 
always ahgned along x. The minimum transverse spin 
variance is 



_ - 

_L,min 2 



1 r 



{Si 



{si)Y + {{s.,Sy}Y 



(17) 



where the expectation values (. . .) represent the average 
over the quantum state and over the random variable D. 
The notation {, } stands for the anticommutator. Intro- 
ducing quantities A and -B, 



A 
B 



^ + A-^ 



52 



(18) 
(19) 

(20) 



To derive the scalings ([T2|) and ([141) ^ ^'^'^ to have a phys- 
ical insight, it is convenient to reason in terms of the 
phases of the operators a and h: 



a — 

h = 



(21) 
(22) 



where Na = (Ya and Nb = b%. This is a legitimate 
representation as long as the condensate modes have a 
negligible probability of being empty By neglect- 

ing the fluctuations of their modulus, the collective spin 
components Sx-, Sy are simply given by 



So, 



N 
Re — e 
2 

Im — e 
2 



i{ea-ei) 



(23) 
(24) 



At t = 0, the phase difference {9a — Ob) has zero mean 
and root mean square fluctuations that scale as 1/ \/N. 
Both Sy and 5*2 scale as VTV- As a consequence, for N 
large we can expand the exponentials in To 
lowest order we then have 



S 



— ^ 



2 ^ 



0, s. 



Na-m 



(25) 



Sy and 5*2 are then simply proportional to the position 
operator Q and the momentum operator P of a fictitious 
free particle. As we will see, the expansions (|25p remain 
valid for times (pgt/h) <C a/a". The squeezing occurs be- 
cause in a given realization of the experiment Sy becomes 
an enlarged copy of 5*2. Indeed after the pulse, for t > 0, 
from the Heisenberg equations of motions for the phase 
operators, with x — 9l{KV)i one has 



%m = (^a-0b)(O+) 



9t_ 

hv 



[2^2 + D] . (26) 



As the squeezing dynamics goes on, Sy that was initially 
of the same order as 5*2, grows linearly in time while 
Sz stays constant. Correspondingly A (x and B oc t. 
Since {{Sy{0), Sz}) = 0, one actually has at all times: 



Sy = Sy{0) + Sl^^''t; 

A 2 

^ = at ; 

- = fit, 



(27) 
(28) 

(29) 



where we have introduced the time independent operator 
and coefficients 



lead 



P9_ 
2h 



'2Sz + D] 
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{{S'^^^Sz}) 
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(30) 
(31) 
(32) 



Using the fact that from 1^5^) S^ ^ N/2 for (pgt/h) < 
VN, and expanding the expression (I^Ul) for t ^ h/{pg), 
we have in the thermodynamic limit 



e(t) = 1 

r— >-oo 



^ + ^l + o(t-) 



Using equations with en 

obtain in the long time limit 



(33) 
1, we finally 



^ (^) — ^noisc 
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pgt 



[l + O (e„oiso)] + O 



\iP9t)' 



(34) 



Best squeezing and close-to-best time 



According to (IMl) . the best squeezing in the thermo- 
dynamic limit to leading order in Cnoiso is 



Sn 



liniTv- 



• N 



(35) 



Remarkably, the best squeezing ([55t only involves the 
part of the phase difference D that is not proportional to 
^2. To understand physically this result we rewrite 



Sn 



^2 ((glead)2^^g2^,^|glcad^^J/2^ 
((5l-d)2)(52) 



(36) 

Due to interactions, through the phase difference (^5]) . 
S^y'^^ is proportional to (25*2 -I- I?). In the absence of the 
term D this allows a perfect cancellation between the 
correlation {Sl'""^Szf and the product ((4°^'^)^)(5'2) in 
([36]) leading to ^^j^ = in the limit A^ — >• oo. In presence 
of D, this is not possible and Cmin ^ non-zero limit. 

Considering the next to leading order in the time ex- 
pansion of the squeezing parameter Eq. p4p . the best 
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squeezing is reached in an infinite time (in the thermo- 
dynamic hmit). However as we will see £,'^{t) is quite 
flat around its minimum, and it suffices to determine a 
"close-to-best" squeezing time t,^ defined as 



Then, according to ([M)) . is given by 



2 

min 



(37) 



(38) 



The close-to-best squeezing time i,, is thus very simply 
related to the best squeezing 4nin- "^^^ important point 
is that trf is finite (non-infinite and non-zero) in the ther- 
modynamic limit. 
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FIG. 3: (Color online) Left: Squeezing parameter as a 
function of time for enoiso = 0.05 (full line) and enoisc = 
(dotted line). Determinant of the covariance matrix M as a 
function of time (dashed line). Simulation with 10* realiza- 
tions. iV = 4 X 10*, X = 0.022577s-^ Right: Isocontours of 
W{y,z) for f = 0, t = 2 X lO^^s and t = 2.5 x lO^^s. 



2. Geometrical interpretation 

We give here a geometrical and pictorial interpretation 
to the squeezing process in the thermodynamic limit. Let 
us introduce the rescaled transverse spin components 



Y 



Si, 



(39) 



that verify (Y) = (Z) = and (Y^) = (Z^) = 1. At short 
times, the Wigner function representing the probability 
distribution of Y and Z is approximately Gaussian 



W{y,z)(xexp -^{y,z)M' 
where M is the covariance matrix 
M = 



{{Y,Z}) {Z') 



(40) 



(41) 



We can represent graphically the fluctuations of Y and 
Z by drawing isocontours of W{y,z). Let us introduce 
the eigenvalues of M 



Ai,2 = l±i({r,Z}) 



(42) 



and a rotated coordinate system y' — z' aligned with the 
eigenvectors of AI: 

y'^iz + y)/V2 ; z' = (z - y)/V2 . (43) 
The points in the y' z' plane such that 

/2 /2 

^ + ^ = 1 (44) 

form an ellipse whose semi-axis gives the mean square 
fluctuations of Y' and Z' that are time dependent linear 
combinations of Sy and Sz. The ellipse surface iSoUipsc — 
7r(AiA2)^/^ divided by tt is equal to the square root of the 



determinant of M and is thus asymptotically equivalent 
to ^inin according to ((55)) : 



= for t» A 

TT pg 



(45) 



In Fig 12] (left) we show the time dependence of the squeez- 
ing parameter in presence and in absence of decoher- 
ence. On the same plot we show the determinant of M 
that asymptotically gives the value of f^j^ ([5^ . In Fig [3] 
(right) we show the isocontours of W{y, z) defined by ([33]) 
at different times for the case with decoherence. As the 
dynamics goes on, Y and Z become more and more corre- 
lated and the ellipse shrinks. In the absence of dephasing 
and in the thermodynamic limit the ellipse would collapse 
into a segment in the y — z direction. In the presence 
of decoherence the process is "blocked" and the ellipse 
keeps a finite width with a limit area iSdiipsc = ""enoiso- 



B. Exact solution of the dephasing model 

The dephasing model (|15p is exactly solvable. One 
first writes the Heisenberg equations of motion for a and 
b, e.g. 



(46) 



Then one uses the fact that Sz is a constant of motion 
to integrate the equations [l^, [ill . One obtains (see also 

My- 



with 



A 
B 
C 



N{N - 1) 



(J2 



l-e-2(xt)^<^'>(cos2xt) 



N-2 



(47) 



(48) 



sinxt (cosxt)^-' (49) 



(5.)^fe-^(>^*)^<^^)(cosxO 



N-l 



(50) 
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A first application of tfie exact solution (|47)) - (|50p is to 
determine the best squeezing ^^j^^ in the thermodynamic 
limit, to all orders in the dephasing parameter enoisc- To 
this aim we take the limit — > cxi in (|T7l) - ([5n)) at fixed 
time t, density p and noise parameter enoisc = (D'^)/N. 
We find that A/N, B/N and C/N have a finite limit and 
that 



(51) 



From this solution one gets the best squeezing and the 
close-to-best squeezing time trj (for ?7 enoisc < 1): 



PQ lim. therm. 



^noisc 

1 



inn 



(1 



(52) 
(53) 
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FIG. 4: Squeezing as a function of time for the dephasing 
model (jlSp as given by the exact solution H47|) . — lO", 
{D^)/N = 0.02. Horizontal dashed lines: given by (|X4| 
and (1 -I- ?7)Cmin with rj = 0.1. Vertical dashed lines from left 
to right: Close-to-best squeezing time f,, (|53p . best squeezing 
time t-mi-a (IA3I) and tL deduced from (IA5I) . 



Note that one can obtain ([55|l and (|38p from ([5^ and 
((53)l by linearizing for small enoisc- 

In FiglHwe show ^^(t) as a function of time for a large 
atom number. The curve is indeed quite flat around the 
best squeezing time tmin- There are two solutions to 
equation ([371): trj < tmin and t'^ > tmin- When A oo, 
t,j is finite and given by (|53p. On the other hand, as 
we show in Appendix [XI t^ain diverges as N^^'^ and tj^ 
diverges as A^/^. Knowing the asymptotic behavior of 
^min, by introducing appropriate rescalings of the time 
variable as in [1^, it is also possible to obtain the first 
finite size correction to Cmin- This is given in equation 
of Appendix [S 



C. Squeezing in the weak dephasing limit 

We can use the exact solution of the dephasing model 
(H71) - ([50]) to investigate the weak dephasing limit that we 
define as the limit where the noise D remains bounded 
for A^ — >■ oo (see the footnote before equation ([T| 



constant A^ 



(54) 



In this case the scaling of the best squeezing time with A^ 
is still given by pgtmin/^ « N^^^ as in the case without 
decoherence pT|) and, using the same rescaling of time 
as in [llj, one has 



III. PARTICLE LOSSES 



In this section we consider particle losses that are an in- 
trinsic source of decoherence in condensed gases. Among 
those, one-body losses are due to collisions of condensate 
atoms with residual hot atoms due to imperfect vacuum. 
More fundamental in dense samples are three-body losses 
where, after a three-body collision, two atoms form a 
molecule and the third atom takes away the energy to 
fulfill energy and momentum conservation. After such 
a collision event the three atoms are lost. Three-body 
losses are present due to the metastable nature of ul- 
tra cold gases, whose real ground state at such low tem- 
peratures would be a solid and whose gaseous phase is 
maintained because the sample is very dilute. Finally 
two-body losses can also be present, caused by two-body 
collisions that change the internal state of the atoms. For 
a trapped gas, we have shown theoretically [l^l that the 
best achievable squeezing within a two-mode model at 
zero temperature in presence of one, two and three-body 
losses can in principle be very large (squeezing parame- 
ter of the order of 10"^) provided that the harmonic 
trapping potential is optimized and a careful choice of 
the internal state of the atoms is made. To realize such 
conditions that minimize losses remains however an ex- 
perimental challenge. 



2 

min 



Pgtmin 



3^/^ 1 

2 A2/3 + N 

3i/6Ari/3„^+o(l) 
4 



(56) 



In this section we recall the main results of [l^ con- 
cerning the squeezing in presence of particle losses, and 
we use these results to show an analogy between the effect 
of the losses and the effect of the dephasing Hamiltonian 
psp in the thermodynamic limit. 
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Monte Carlo wave functions 



t-t 



We consider two spatially separated, symmetric con- 
densates. For a more general treatment, please refer to 
[ll|- Initially the system is in the eigenstate of with 
maximal eigenvalue N/2. Besides the non-linear Hamil- 
tonian for the two bosonic modes a and b given by 



Hni = hxSl with {dN^^J-a)^ 



(57) 



we include one, two and three-body losses. Due to the 
losses, the system is "open" and we shall describe it with 
a density operator that obeys a Master Equation of the 
Lindblad form [To| . In the interaction picture with re- 
spect to Hnl- 



Jm rn ~ 



rp} 



where Cf — a,b for e — a,b and: 

b = 



(58) 



(59) 



Here the operators a and b are in the Schrodinger picture. 
The m-body loss rates j'-"^'' are defined in terms of the 
so-called rate constants Km as 
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(m) _ 



dV|</)(r)| 



2m 



(60) 



where ^(r) is the condensate wave function in mode a or 
b for the initial atom number (weak loss approximation) , 
so that for example: 



dt 



Ki + K2Na / dV|0(r)|4-f 



KsN! / dV|0(r)|6 Na 



(61) 



It is convenient to rephrase the Master Equation (1551) in 
terms of Monte Carlo wave functions [17]] . In this picture 
pure states evolve deterministically under the influence 
of an effective Hamiltonian iJcff acting during time inter- 
vals Ti — ti ~ ti^i separated by random quantum jumps 
(described by the jump operators S^) occurring at times 
ti as illustrated in FigIS] 

ih 



5, = ^7^4". (62) 



e—a,b 



More precisely the evolution of the non-normalized state 
vector {ipit)) in between quantum jumps is given by 



ih-\m) = Hes{t)\m) 



(63) 



and the effect of a quantum jump in the component at 
time ti is 



(64) 



t 



k-l 



FIG. 5: Time sequence of deterministic evolution periods and 
quantum jumps for a single Monte Carlo wave function \ip). 



Quantum averages of any atomic observables O are 
obtained by summing over all the possible trajectories of 
the non- normalized state vector pjj] : 

(0)-E/ rfiidi2---dife ^ (^(i)|0|i/;(t)). (65) 



B. Losses randomly kick the relative phase 

Let us consider the action of a quantum jump over a 
phase state ^ with N atoms, for example the loss of one 
particle in state a or 6 at time t: 



CaitMN = \/ye-*4*e'*|</.-xt/2)Ar_i 



(66) 



Cb{tMM = ^ye-*^*e-^^|</. + xt/2)iv-i. (67) 

Under the action of a jump a phase state remains a phase 
state. On the other hand, the relative phase is shifted 
by a random amount that depends on the time of the 
jump and has a random sign depending on whether the 
jump was in a or b. This behavior is illustrated in FigjS] 
taken from ilOj where we plot the modulus squared of a 
relative phase distribution amplitude c((/), t) &t t = Q and 
t ~ 27r/x for three single Monte Carlo realizations. At 
this particular time (second revival time) the coherent 
evolution due to i?ni has no effect and we can isolate the 
action of the losses. In the case of FiglBl as xt is of the 
order of unity, the shift in the relative phase due a single 
jump is large. In the case of squeezing with iV 3> 1, 
Xt ^ 1, the shifts due to single quantum jumps are on 
the contrary very small. These shifts nevertheless limit 
the maximum squeezing achievable as we shall see. 



C. Spin squeezing limit and the lost fraction 

Let us consider the case of one-body losses only with 
a loss rate constant 7 equal in states a and b. In this 
case the effective Hamiltonian does not depend on time; 
referring to the time sequence in FigjSl we have explicitly 

\^{t)) = e-'^-(*-*'=)/'^5,,(tfe)e-'^^--'=/''.S,,_,(tfe_i) . . . 

...S,,{h)e-^"'-^-'^\m), (68) 

and there is an explicit analytical solution for the gener- 
ated spin squeezing as a function of time [ll|, • Here 
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In Fig 17] we show the squeezin g in presence of losses and 
we compare the exact solution |12| with the approximate 
expression ([72| valid at long times and in thermodynamic 
limit. 



>u- 10 




FIG. 6: Relative phase probability distribution at time t — 
and at time t = 2n/x in three single realizations Monte Carlo 
realization. From upper left to lower right the wave function 
has experienced 0, 3, 1 and quantum jumps respectively. 
Three-body losses corresponding to = 301 *^Rb atoms in 
F = l,mf — —1 state in separated identical harmonic traps 
with to = 500Hz. Figure taken from p^ . 



FIG. 7: (Color online) Spin squeezing as a function of time 
in presence of one-body losses: 7 — 0.005pg/h, N = 10^. 
Solid black line: Exact analytical solution. Dashed red line: 
Approximate solution (|72l) . Horizontal and vertical dashed 
lines: predictions (|73|) and (|74p respectively. 



we use this solution to find the best squeezing in pres- 
ence of losses in the thermodynamic limit with N ^ 00, 
p and 7 constant. For simplicity we also assume that the 
fraction of lost particles at the relevant time t remains 
small in the thermodynamic limit: 



it = eio 



(69) 



We proceed similarly as we did to derive Ea. ((CT|) . taking 
the thermodynamic limit in the exact solution in presence 
of losses. Restricting for simplicity to the leading order 
in 7i and to the case pgt/h ^ 1, we obtain 



A 



B 











-It) 


hi 


-) 







[1 + 0(7i)] 



(70) 
(71) 
(72) 



where A and B are defined in ([T8l)-(|T9l). Equation (|72|) 
shows that for long times, the squeezing parameter is 
asymptotically equivalent to one third of the lost fraction 
of atoms. Minimizing (|72p with respect to time we obtain 



4 h-y 
^ pg 



2/3 



j_ '-min 

n 



4 h-y 



1/3 



(73) 
(74) 



D. Analogy between dephasing and losses 

Our goal here is to make an analogy between the de- 
phasing model of section |lT] and the present model with 
losses. Indeed the relative phase is perturbed by the 
losses. In presence of one-body losses, a non normalized 
Monte Carlo wave function after k jumps has the form: 



\m) 



II' 



-7U g-7(JV-fc)t^fc 



V = 



N{N ~1)...{N -k + 1) 



1 



(75) 

(76) 
(77) 



1=1 



The factor Af given by ((76)) is the norm squared of the 
wave function that is needed to calculate quantum aver- 
ages, the first line in ([76l) is due to the effective Hamilto- 
nian evolution, (j) is an irrelevant phase and V given by 
(|77p is a random perturbation of the relative phase 20 
that plays the role of the quantity D in Eq. (P5| of the 
dephasing model of section |lTl as appears from the fact 
that exp{—ixtDSz)\(l))N = \4> — XtD/2)jq. Contrarily to 
Z?, V given by ((77)) is time dependent. As detailed in 
Appendix iBl using (I55|) we can calculate (I?^). To first 
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order in we obtain: 



N 



3 



(78) 



We have thus shown that ^^(i) is asymptotically equiv- 
alent to {'D'^)/N as it the case in the dephasing model. 
We can then establish an analogy between the model with 
losses and the dephasing model as summarized in the first 
two columns of the Table HI In Figl8]we show a compari- 
son of {T>'^)/N obtained from a Monte Carlo simulation, 
from the exact expression (jB7[) . and from the approxi- 
mate expression ((78|) valid to first order in eioss = it- 



0.06 



0.04- 



A 
V 



0.02 




FIG. 8: (Color online) Main plot: {■D^}/N given by JZZl as a 
function of time for 7 — 0.005pg/h, N — 10^. Black circles: 
Numerical simulation with 5 x 10^ realizations. Dashed blue 
line: Approximate expression (|78() . Solid blue line: Exact 
expression (|B7|) . Inset: Squeezing as a function of time. Solid 
red line: Exact analytical solution. Black circles: Numerical 
simulation with 5 x 10^ realizations. 



number. This optimum squeezing in a trap in presence of 
one, two and three-body losses has a simple expression as 
a function of the s-wave scattering length Uaa = afeh ~ a 
and the rate constants Km, tti = 1,2,3: 



2 ^ /5V3M 
^°P' I 287r ha 



2/3 



(K1K3) + K2 



2/3 



(80) 



In Figiniwe show the spin squeezing minimized over time 
^^ijj as a function of N in presence of one, two and three- 
body losses. The trap frequency is optimized for each 
value of iV fil: 



219/1275/12^5/6 ^ ^1/2 



^opt 



151/3 



iVl/3 



K3 



5/12 



(81) 



We note that the minimum squeezing, instead of going 
to zero as in two-mode model without decoherence (red 
line) tends to a finite non-zero value for TV — > cx) given 
by (gOl). 
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E. Optimum squeezing in a harmonic trap 

A similar analysis can be performed for a trapped sys- 
tem, where we now consider the two components a and 
b in identical and spatially separated harmonic traps. In 
particular equation (19) of |12| is very similar to ([7^ 
except that in (19) of ^1^ we have Nx instead of pg in 
([7^ . Another difference is that the more general equa- 
tion (19) in [l^l that includes two and three-body losses 
besides one-body losses, is derived performing an approx- 
imation on the effective Hamiltonian: the constant loss 
rate approximation (25| 

= f T^-^ct-C ^ - E f T^-^^r • (79) 

e— a,fe e—a^b 

By using (19) of [l3| and the Thomas-Fermi profiles 
of the condensate wave functions, one can optimize the 
squeezing with respect to time, trap frequency and atom 



FIG. 9: (Color online) Blue line: Best squeezing 5,^;^ as 
a function of A'^ in presence of one, two and three-body 
losses, obtained in hj in the constant loss rate approxi- 
mation. The trap frequency is optimized for each value of 
N according to (|8ip . Dashed blue line: Optimum squeez- 
ing for large N given by the analytical prediction (|80|l . Red 
line: Best squeezing without losses for comparison. Param- 
eters: a = 5.32nm, if^^' = 0.1s-\ K^^'> = 2 x 10-^^m^s-\ 
K^^^ = 18 X 10"^'^m^'s~\ This figure is taken from 



IV. FINITE TEMPERATURE 

The multimode nature of the atomic field and the 
population in the excited modes at non-zero tempera- 
ture have important consequences on the squeezing and 
change the scaling laws with respect to the two-mode 
case (fTTI) . In [1^ we use a powerful formulation of the 
Bogoliubov theory in terms of the time dependent con- 
densate phase operator [2l[ to perform a multimode 
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treatment of the squeezing generation in a condensed gas 
in the homogeneous case. We show that the best squeez- 
ing ^min has a finite non-zero value in the thermodynamic 
hmit and we calculate this value analytically. 



A. Multimode description 

We consider a discretized model on a lattice with unit 
cell of volume dV, within a volume V with periodic 
boundary conditions [2l[. The Hamiltonian after the 
pulse for component a (and similarly for b) reads: 



2m 



The fields have commutators 



dV 



(82) 



(83) 



with 1.1,1/ — a or 6, and ak(6k) is the amplitude of ^a.b 
over the plane wave of momentum k. We assume identi- 
cal interactions in states a and b with a coupling constant 
g — ATrfi^a/m where a is the s-wave scattering length 
while states a and b do not interact {gab — 0). Note that 
the coupling constant in Ha.b should actually be a bare 
coupling constant gQ different from the effective coupling 
constant g, but this difference can be made small in the 
present weakly interacting regime (pa^ Y^^ <C 1 by choos- 
ing a lattice spacing much larger than a but still much 
smaller than the healing length cx l/^/pd [22]. 

In terms of the fields, the collective spin components 
are 



S,+iSy = dVj:,tPl{r)Mr) 



(84) 



with = dy^ri/;t(r)V'^(r), = a,b. 

For low or high values of ksT/ pg we find the asymp- 
totic behaviors 



f 2 (T=0) 

kBT<^pg: ^02= ~ ^Siis^ = 0.0234 . 



ksT pg: C 




(85) 



(86) 



B. Best squeezing and close-to-best time 

Performing a double expansion for large N and small 
non-condensed fraction eboe, 



CBog = — ^ < 1 



(87) 



we find that to first order the system effectively behaves 
as in the dephasing model presented in section |TT1 Indeed 
the component Sy of the spin develops a term that is 
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FIG. 10: (Color online) (solid line) and before-pulse non- 
condensed fraction {Nnc)/N (dashed line), both divided by 
^fpa?, as functions of ksT/pg. Both quantities are obtained 
analytically within a Bogoliubov framework, see H92[) for ^min- 



proportional to the condensate relative phase 9a — df,, 
and the relative phase evolves as 

gt 



% -Ob= {6a - 0fc)(O+) - ^ [{Na - Nb) 



The dephasing parameter D^h in 
population in the excited modes: 



Ah] . (88) 
is related to the 



D 



th 



where n^k — c^a^Cak and nf,k = c^k'^bk are the occupation 
number operators of the Bogoliubov modes. Note that 
these modes are in a non-equilibrium state since the zero 
relative phase state between the condensates in a and b 
is prepared at i = 0"*" by applying a sudden 7r/2 pulse to 
the gas (that was initially at thermal equilibrium in state 
a). Uk and Vk are the usual Bogoliubov functions 



{Uk + Vkf {UaV 



nbk 



(89) 



Uk + Vk 



Ek 



Ek + P9 



1/4 



h^k^ 

2m 



(90) 



p/2 = N/2V being the spatial density in each single com- 
ponent a or 6 after the pulse. Dth ffuctuates from one 
realization to the other because Uak and n{,k depend (in 
Heisenberg picture) on the creation and annihilation op- 
erators of the Bogoliubov modes before the pulse, and 
the initial state of the gas has thermal fluctuations. 

To first order in N and in enog, the best squeezing 
parameter and the close-to-best squeezing time are given 
by H 



(A\> 



pgtv 



1 



N ' h 

An explicit calculation [l^ gives: 

Smin 



(91) 
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Particle Losses 


Dephasing model 


Multimode T 7^ 




{Oa - 0b){t) = {Oa ~ eb){0+) - xt [25. + D] 


(Sa - = {Oa - eb){Q+) - Xt [25. + Ah] 


T> from quantum jumps (|77|) 


D from a dephasing Hamiltonian (|15|) 


Dth from excited modes population (|89|l 


en) - ^""'^ 

pgt/h>l N 


^^(^) ~ ^""'^ 

P9t/h>i TV 


P9t/h>l 


3 3 


— tnoise 


^ = y^F{kBT/pg) ^ ^Bog 

iV kBT>pg 



TABLE I: Correspondence table among the different physical models. The results are valid in the thermodynamic limit and to 
first order in the small parameters eioss (f5^ . enoise (fT5)l and esog ((57)) . 



where now n^,"^ are Bose mean occupation numbers of 
Bogoliubov modes in state a before the pulse 



ek = ^Ek{Ek + 2pg). (93) 



Sk = t/fe + Vfe and si. 



(0) _ ^(0) 



(0) 



where t/^*^' and Vj, 



(0) 



are the Bogoliubov functions in internal state a before 
the pulse 



(0) 



Ek + 2pg 



1/4 



2 7,2 



Ek 



2m 



(94) 



A consequence of (1921) is that the squeezing divided by 
\fpa^ is an universal function of kgT I pg- 

C. Spin squeezing and non-condensed fraction 

Within our treatment, valid for large TV and T ^ Tc, 
we find that the best squeezing ^^j^^ (|92l) is always lower 
than the before-pulse non-condensed fraction. This is 
shown in FigfTOl 

Finally in Table IIV Al we can complete the correspon- 
dence table between dephasing noise, losses and non-zero 
temperature effects. 



V. A FEW WORDS ABOUT EXPERIMENTS 

Two recent experiments demonstrated spin squeez- 
ing in Bose-condensed bimodal condensates of rubidium 
atoms in internal states la) = \F = \^mp = =pl) and 
\b) = \F = 2, nip = ±1) Due to the fact that the 

three scattering lengths characterizing the interactions 
between atoms are very close 



ffaa — gab — Qbh 



the effective two-mode nonlinearity x in H^i (HJ is very 
small when the condensates a and b overlap spatially. 
This gives the possibility to tune the non-linearity [11| 
by either controlling the spatial overlap between the two 
species as done [6j or by using a Feshbach resonance 
changing the inter-species coupling constant gab as done 
in 4]. A fundamental source of decoherence in these 
experiments is two-body losses in the state F = 2. A 
theoretical analysis [ll| shows that in typical experimen- 
tal conditions these losses limit the squeezing to about 
~ 6 X 10~^. For cold samples this limit is above the 
limit imposed by non-zero temperature. In the case of @ 
we could explain in detail the squeezing results using a 
zero-temperature model including spatial dynamics and 
including particle losses and technical noise (dephasing 
noise) as sources of decoherence, the latter being domi- 
nant [6, 18]. 

The perspective of studying experimentally the scaling 
of squeezing in a controlled decoherence environment, e.g. 
preparing the sample at different temperatures, is fasci- 
nating and challenging. 



VI. CONCLUSIONS 

We have considered a scheme to create spin squeez- 
ing using interactions in Bose-condensed gas with two 
internal states 0, 0, [1| ■ The squeezing is created dynam- 
ically after a 7r/2 pulse applied on the system initially 
at equilibrium in one internal state. We have reviewed 
the ultimate limits of this squeezing scheme imposed by 
particle losses and non-zero temperature based on our re- 
cent works [ll] and [ll] and we have extracted a simple 
physical picture of how decoherence acts in the system. 
An important result is that contrarily to the case without 
decoherence the squeezing parameter minimized over 
(95) time ^^j^^ has a finite non-zero value in the thermody- 
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namic limit, that we determine analytically. Finally we 
have shown that the physics of spin squeezing in presence 
of losses or at non-zero temperature can be caught by a 
simple dephasing model (also considered in [l^) that we 
have solved exactly and studied in details. 



Appendix A: Times t,nin and t'^ in the dephasing 
model 



The exact solution (|T7)) -([5(I 1) allows to determine how 
the best squeezing time tmin diverges in the thermody- 
namic limit. We first found numerically that it diverges 
as N^^'^. We then introduce the rescaled time 9 such that 



P9 



(Al) 



and 



Expanding the functions cos (2xt) 
exp[-2(xi)^(£)^)] up to terms Oll/N) included in 

A] linearizing sin(xt) and expanding cos^~^(xi) and 
exp[-(xi)^(-D^)/2] up to terms 0{l/N'^/'^) included in 

B] and expanding cos'^~^{xt) and exp[— (xi)^(-D'^)/2] 
up to terms 0(1/A^^/^) included in C, we obtain 



Appendix B: Calculation of (Vi^) in the lossy model 

From the definition of T) ([77)1 and the expression of 
a quantum average (j65|) in the Monte Carlo wavefunc- 
tion method, using the expression of the norm squared 
of \tp{t)) of (I76|, we obtain 



h E {\tnA\i{e-A. (Bl) 

»7i,...,i)fc = ±l \ 1=1 / \i=l / 

We have introduced the random variables rji = +1 for 
ei = b and iji = —1 for e; = a, and we have used the 
fact that the integrand is a symmetric function of the 
jump times ti to extend time integration from the ordered 
domain < ti < . . . < tfc to the hypercube [0,^]*^ (also 
dividing by fc!). The notation (^) represents the usual 
binomial coefficient A^!/[fc!(iV — fc)!]. We first sum over 
the variables rji: 



iVl/2 



noise 
fl2 , 



1 



^^(l-f e„oisc)3 



By minimizing (|A2I) over 9 one obtains 



Pff , lim. therm. 

n 

J.9 lim. therm. 



N 



=(1 



1/4 



l + En 



1/2 



iVl/2 



(1 



3/2 



o 



N 



(A2) 



(A3) 



A4) 



With a similar technique we can determine the diver- 
gence of t'^ that is the second solution of equation ([37)) 
t'^ > imin- When TV -> cx), diverge as iV^/^. To cal- 
culate the prefactor we introduce again a rescaled time 
pgt/h = 0'[N/{1 + e„oisc)]^/^ and take the large N limit 
in ([TT)) to obtain 



1 - 



sinh 9 



a 



(A5) 



Solving the transcendental equation (1 + ?/)C'^ = Cmin o^^^ 
finds the large A'' approximation to t'^ . For 77 ^ 1 and no 
constraint on the ratio 77/enoisc, we obtain 



9'2 



2f] 



2>? 
3e„oi3 



(A6) 



4 E iT.^A =E^' 



2'= 

Jji ,...,Jjfc=±l \i=l / 1=1 

then we perform the temporal integration to obtain 

N 



(B2) 



E 
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,fe^2 



k J Iq 



k=l 

* 1 - e-^* 



/o = / dii e-^*^ = 
Jo 

h = f dhtle--"' 
Jo 



(B3) 
(B4) 



= ;^{2-e-^*[2 + 27i+(7t)2]} , (B5) 



(B6) 



Taking the derivative with respect to u of the binomial 

• k > 



identity J2i=Q i^)^'' = (1 + '")^: we get the final expres- 



sion 



N ■ 

Expanding I2 for small 7t gives as expected 
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